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| Classical mechanics (CM), like quantum mechanics (QM), can have an operatorial 
formulation. This was pioneered by Koopman and von Neumann (KvN) in the 30 's. 
They basically formalized, via the introduction of a classical Hilbert space, earlier work 
of Liouville who had shown that the classical time evolution can take place via an 
operator, nowadays known as the Liouville operator. In this paper we study how to 



a: 



perform the coupling of a point particle to a gauge field in the KvN version of CM. 
So we basically implement at the classical operatorial level the analog of the minimal 
coupling of QM. We show that, differently than in QM, not only the momenta but 
also other variables have to be coupled to the gauge field. We also analyze in details 
how the gauge invariance manifests itself in the Hilbert space of KvN and indicate 
the differences with QM. As an application of the KvN method we study the Landau 
problem proving that there are many more degeneracies at the classical operatorial 
level than at the quantum one. As a second example we go through the Aharonov- 
Bohm phenomenon showing that, at the quantum level, this phenomenon manifests its 
effects on the spectrum of the quantum Hamiltonian while at the classical level there 
is no effect whatsoever on the spectrum of the Liouville operator. 
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1 Introduction 



It is well-known that in classical statistical mechanics the evolution of probability 
densities p(q, p) is given by the Liouville equation 

dtp(q,p) = -Lp(q,p) (1.1) 

where L = d p Hd q — d q Hd p and H(q,p) is the Hamiltonian of the system. The p{q,p) 
must only be integrable functions (i.e. p G L 1 -space) because of their meaning as 
probability densities {i.e. J p dqdp < oo). As they are only L 1 -functions they do not 
make up an Hilbert space. In 1931 KvN [|I|] postulated the same evolution equation for 
complex distributions if>(q,p) making up an L 2 Hilbert space: 

d t i[)(q,p) = -Lif)(q,p) (1.2) 

If we postulate eq. (|1.2|) for if)(q,p), then it is easy to prove that functions p of the 
form: 

P=H 2 (1.3) 

evolves with the same equation as if>. This is so because the operator L contains only 
first order derivatives. This is not what happens in QM where the evolution of the 
if>(q) is via the Schrodinger operator H while the one of the associated p = \ip\ 2 is via 
a totally different operator. The reason is that the Schrodinger operator H, differently 
than the Liouville operator L, contains second order derivatives. By postulating the 
relations (O) and ( |1.2| ) for the ip, KvN managed to build an operatorial formulation 



for classical mechanics (CM) equipped with an Hilbert space structure and producing 
the same results as the Liouville formulation. Of course there are phases in a complex 
if) and one could wonder which is their role in CM. This problem is addressed in ref.@. 
We will briefly review the KvN formalism in section 2 of this paper. 

The question we want to address in this paper is how to couple a gauge field to the 
point particle degrees of freedom once we work in the operatorial formulation of CM. 
We know that in QM there is a simple rule known as "minimal coupling" (MC) 0: 

g _ 

p — > p A (1.4) 

c 

-» 6 -» 

where A is the gauge field. This rule says that it is enough to replace p with p A in 

c 

g 

the Hamiltonian H(q,p) — ► H(q,p — A) in order to get the interaction of the particle 

c 

H d 

with the gauge field and then represent p operatorialy as — — . Since the operatorial 

i ox _ 

formulation of CM has an Hilbert space structure like QM and it has operators like L 
analog of the Schrodinger H operator of QM, we would like to find out the "minimal 
coupling" rules which would transform the L without a gauge field interaction into the 
one L A with interaction. 
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Let us start by just using the rule (|1.4j) inside the H which appears in L: 

L = (d p H)d q — (d q H)dp (1.5) 
d 

without modifying the derivative — . We can work out the simple case of a particle 

dp 

moving under a constant magnetic field directed along z. Let us choose the gauge field 



as: 



Using the MC 0: 



the Hamiltonian becomes 



A x = 

A y = Bx (1.6) 
A z = 



Py — ► Py ~ ~ x (1-7) 



Pi 1 ( eB V Pi c 



2m 2m \ c J 2m 
and L would turn into 

1 = ?ljL + L(p eB x]( 9 I eB 9 ) I Pz 9 (19) 
A mdx m \ v c )\dy c dp x ) mdz 

If we compare this Liouville operator with the one containing no interaction with the 
magnetic field which is 

t = r.B +bL a + p.a 

m ox moy m oz 
we see that the tricks to pass from ( |1.10|) to ( |1.9|) are the substitutions: 

eB 

Py > Py ~ x 

d d eB d 

+ 



dy dy c dp. 



X 



These are the MC rules for the Liouville operator in the case of a constant magnetic 
field. 

For the Schrodinger operator the MC rules would have been the quantum operatorial 
version of just the first relation of (|1 . 1 1| ) 

hd_^hd__e_B_ x 

i dy i dy c 

h d 

which would turn into the second one if we had represented x = — — — . Somehow in 

i dp x 

QM the two rules (|1 . 1 1|) would be just one and the same differently than in CM where 
we cannot identify p y with — . 



In section 3 of this paper we will generalize the MC rule ( 1.11 ) to the case of an 
arbitrary magnetic field. Those rules, even if derived from the same simple principle as 
above, will involve various complicated combinations of the derivative operators. We 
shall show anyhow that those complicated combinations could be put in a very simple 
and illuminating form using the concept of superfield which naturally appears in a 
functional approach 0] to the KvN theory. This functional approach will be briefly 
reviewed in section 2. 

Let us now go back to our derivation of the minimal coupling rules for L. If in the 
original H there were a potential V(y) then the L of eq. (|1.5f) would have contained 
a derivative with respect to p y . Now in (|1.7|) we changed p y to get the MC and so 

d d 
one would be led to conclude that also — — has to be changed in L and not just — — 

dp y dy 



as it appears in (|1 . 1 1| ) . Actually this is not the case . We will see in section 3 that 
the derivatives with respect to the momenta never have to be changed in the MC 
correspondence. 

The reader, realizing that we have an Hilbert space, may wonder if the gauge invari- 
ance manifests itself on the states via a phase like it does in QM. This is so only in a 
particular representation. The issue of gauge invariance will be throughly examined in 
section 3 of this paper. 

In section 4 we will apply the KvN theory to the Landau problem. We will compare 
the results to the quantum case and show that there are many more degeneracies at 
the classical level than at the quantum one. Finally in section 5 we will study the 
Aharonov-Bohm phenomenon both at the CM level via the KvN formalism and at 
the QM level via the Schrodinger equation. We will show that at the quantum level 
the spectrum of the Schrodinger Hamiltonian is changed by the presence of the gauge 
potential while the spectrum of the classical Liouville operator is left unchanged. This 
spectrum drives the motion of the system at the classical level as it is explained in 
ref. [0]. Further calculational details omitted in the various sections of this paper are 
confined to few appendices which conclude the paper. 



2 Functional Approach to the KvN Theory 

As we indicated in the previous section, KvN postulated for the Hilbert space states 
ip of CM the same evolution ( |1.2| ) as for the probability densities p (|1 . 1|) . So the 
propagation kernel for ip will be the same as the propagation kernel for p. This last 
one has an immediate physical meaning being the transition probability P((p a t\(pQto) 
of finding the particle in the phase space point (p a = (q,p) at time t if it was in 
configuration (p^ at time t . In CM this P((p a t\(fQto) is nothing else than a Dirac delta 

P(^% %) = % a -^(t;^o)] (2-1) 
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where is the classical solution of the Hamilton equations of motion up a = uj ab -^—^ 

(u ab symplectic matrix) with initial condition (ipo,to). We know that in general, if 
we have a probability P(f\i) to go from configuration (z) to configuration (/), the 
following decomposition holds: 

P(f\i) =Y, P (f\k N . 1 )P(k N „ 1 \k N _ 2 ) ■ • -P^ll) (2.2) 
In our case this becomes 

N 

P(^|^o) = Km n / d 2n ^ 2n) {v] ~ mM-ih-i)} (2.3) 

where we have sliced the interval of time t — to in N intervals labelled by tj. In the 
continuum limit we could formally write the relation above as 

P{<pt\<pfa) = Jv<p - <fc(t)] (2.4) 

where T>(p is a functional integration and 5[ ] a functional Dirac delta. As the 0^(t) 
are the solutions of the Hamilton equation, we could rewrite the S[ } in (2.4) as 



5[cp a (t) - (f> a d (t)] = 6[<p a - u ab d b H}det[5 a b d t - d b (oo ad d d H)} (2.5) 



Next we could Fourier transform the Dirac delta on the RHS of ( p.5|) introducing 2n 
extra variables A a , and we could exponentiate the determinant on the RHS of ( p.5|) 
using An anticommuting variables c a , c a . The final result is the following 



P{^ a t\^to) = I V"(pVXDcVcexp\i [ dtC 

where 



(2.6) 



C = \ a ip a + %c a c a - X a oo ab d b H - tc a u ad (d d d b H)c b (2.7) 

and with T>" we indicate that the integration is over paths with fixed end points in (p. 
All this is described in many more details in ref. H. Associated to this L there is an 
Hamiltonian which is 

H = X a uJ ab d b H + ic a u ad (d d d b H)c b (2.8) 

We notice that, instead of just the original 2n phase space coordinates (p a , we now have 
8n variables (ip a ,X a ,c a ,c a ) whose geometrical meaning has been studied in details in 
ref. ||^ We will indicate with M. the original phase space coordinatized by ip a and 
with M. the space coordinatized by ((p a , \ a ,c a ,c a ). What is nice is that all these 8n 
variables can be put together in a single object known as superfield in the literature 
on supersymmetry. In order to construct it we first enlarge the standard time t to two 
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Grassmannian partners 9, 9 which make, together with the time t, what is known as 
superspace and then we build the following object 

$ a (t, 9, 9) = V a (t) + 9c a (t) + 9u ab c b (t) + i99u ab X b (t) (2.9) 

It is easy to prove that 



. J d9d9H[$]=H (2.10) 



where H is the usual Hamiltonian of classical mechanics in M. but where we have 
replaced the standard bosonic variables (p a with the superfield variables $ a . By ex- 
panding H(Q) in 9 and 9 it is straightforward to prove the following formula 

H($ a ) = H((p a ) + 9N H -9N H + i99H (2.11) 

where the precise form of N H , N H is not necessary in this paper and can be found in 
ref. [Q. The same steps we did above for the Hamiltonian can be done for any function 
0(ip) of the phase space A4, i.e. replace ip with <3> and expand 0($) in 9, 9 

0($ a ) =0(ip a ) + 9N o -9N o + i99O (2.12) 

where O = X a uj ab d b O + ic a u ad {d ( idi ) 0)d > . The index "a" in (p a indicates either the first 
n configurational variables q or the second n momentum variables p, so in the case of 
n = 1 the substitution in H on the LHS of ( |2.10D is: 

q > $1 = q + 9c q + 9c p + i99\ p (2.13) 

p > $P = p + Q C P _ fog - l99\ q (2.14) 

We have put the index q (or p) on c, c, A just to indicate that we refer to the first n 
(or the second n) of the c, c, A variables. A relation analog to (|2.10| ) holds also at the 
Lagrangian level 

i J d9d9L[$] = C + (s.t.) (2.15) 
where (s.t.) is a surface term of the form 

(s.t.) = --(\ p p + ie p c?) (2.16) 

In this paper we will be interested only in the non- Grassmannian set of variables 
among the 8n-ones ((p a , X a , c a , c a ), that means only in (p a , X a . So the C and the H of 
formulas (|2~7|) and ( p.8|) will be reduced tcQ: 

£ B = X a if a - X a uJ ab d b H (2.17) 
^ s = X a u ab d b H (2.18) 



3 The subscript B is for Bosonic but we will neglect it from now on. 
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and the superfields ( |2.13|) and ( |2.14| ) to: 



= q + iOe\ p (2.19) 

<$>p = P -m\ q (2.20) 

KvN postulated for the wave functions ip the same evolution as for the p, so the kernel 
of propagation for the ip: K ((p a t\(pQto) , also known as transition amplitude, will have 
the same path integral expression as the transition probability P(ip a t\ipQt Q ). This last 
has the expression ( |2.tj ) where actually the integration over the c a , c a can be dropped 
because we are not propagating these variables. Their overall integration would give 
one as proved in ref. [[§]. So the equality of P(y9 a t|y9gto) and K(ip a t\(pQto) woulcf] give: 



P(ip a t\^t ) = K(ip a t\^t ) = J WipVX exp iJdtC 



(2.21) 



The two path integrals for P( | ) and for K( | ) are the same if we want to propagate 
probability densities p or wave functions ip both living in the y?-space only. If instead we 
want to propagate the analogous objects living in the (<p, c)-space, then the evolution 
of p(f, c) would be via the £ of eq. ( |2.7D but the evolution of the ip(<p, c) would be via 
a different Lagrangian studied in ref. 0. 

In a manner similar to what is done in field theory we could also introduce the 
generating functional which would have the form 



Z CM [j] = J VipVXexp i J dt [C B + + 



(2.22) 



where j v , j\ are currents associated to the two variables ip, X. We shall now derive the 
operatorial formalism associated to this path integral and compare it with the KvN 
operatorial version of CM given by eq. ( |1.2|) . Having a path integral we can introduce 
the concept of commutator as Feynman did in the quantum case: given two functions 
Oi(ip, X) and 2 (</2, A), let us evaluate the following quantity under our path integral: 

([0 1 ,0 2 ))=^(0 1 (t + e)0 2 (t)-0 2 (t + e)0 1 (t)) (2.23) 

which is defined as the expectation value of the commutator. In our case what we get 
is 

<[<//,</]> = (2.24) 

(& a ,X b ]) = i5 a b (2.25) 

The first commutators confirm that we are doing CM and not QM. In fact they tell 
us that the q commute not only among themselves but also with the p. The second 

4 The reader may be puzzled that the same kernel propagates both tp and tp\ 2 . This is not a 
problem and does not lead to any contradiction as proved in appendix A of ref. [ 2| . 
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commutators instead tell us that the A a are something like the momenta conjugate to 
ip a . In order to satisfy the ( |2.25|) they can be realized as 



A„ = —i 



d 



(2.26) 



Via this operatorial realization of A a also Tt of eq. (|2.18|) can be turned into an operator 

H B — ► H = -iu ab d b Hd a (2.27) 

and it is easy to see that Ti = —iL where L is the Liouville operator of eq. ( |1.5| ). So 
this confirms that the operatorial formalism generated by our path integral is nothing 
else than the KvN one. 

If we specify things to n = 1, the relations (|2.25|) are just 2: 



\p, 



(2.28) 
(2.29) 



and the representation we used for X q in (|2.26| ) is not the only possible one. In fact in 
( f2.26j ) we realized tp a as a multiplicative operator and A a as a derivative one but, using 
the fact that fl2~28"D and (|2T29| ) are two distinct relations, we could have realized q and 
A p as multiplicative operators and A g and p as derivative ones: 



Q 
A, 



V 



—% 



d_ 

dq 



(2.30) 



. d 



In the representation ( ^.26[ ) we would diagonalize (p a 

a \q,p) = v a \<i,p) 

and the kernel ( p.21|) could be formally represented as 

K{(f a t\(f^t ) = (qpt\q p t ) 



(2.31) 



(2.32) 



In the mixed representation ( [2.30|) we would diagonalize instead these other operators: 



q\q, A p ) = q\q, X p ) 



(2.33) 



In this representation we have to do a "partial Fourier transform" in order to go from 
the states \qp) to the \q\ p ) ones . We say "partial Fourier transform" because somehow 
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we are just "replacing" p in the state \qp) with A p which is somehow the momentum 
conjugate to p as can be easily seen from the Lagrangian C B of eq. (|2.17|) . Using the 
transformation formula^ 



{q'p'\q, \ p ) 



'Ik 



S(q — q')exp[ip' X % 



(2.34) 



we can transform ip from the \q,p) basis to the \q, X p ) one: 

1 



dqdp (qX p \q'p) ■ (qp 



i/j(q,p)e- iXpP dp (2.35) 



In the same manner we can transform the kernel (|2.32|) and get: 

K(qX p t\q X po t ) = (qX p t\q X po t ) = 

dq' dp' dq'dp'(qX p t \ q'p't) {q'p't \ q' p' t ) (q' p' t 1 q X po t ) = 



dq' dp' dq' dp 



- iX pP'b( n - a' 



2tt 



5(q-q')K(q'p't'\q' p' t 



271 



<%o - Qo) 



— / dp' dp'exp[-iXpp'}exp[iXp p' }K{qp't\q Q p' t () ) 

Z7T ' 



(2.36) 



This mixed representation is the most useful one in analyzing the issue of how to go from 
this formulation of CM to QM 0. It is also the representation where gauge invariance 
manifests itself via the multiplication by a phase like in QM. We will examine all this 
in details in the next section. 



3 Minimal Coupling in Tt and Gauge Invariance 



In the previous section we have introduced the concept of superfield, formulas Q2.9D - 
( p. 13j )-( [2~T4T )-( |2. 19| )-( |2.20| ) . Using them it is then easy to put in a compact form the 
minimal coupling rules fll.ll ) used in the case of a constant magnetic field. First, via 
the representation ( 2.26|) , the relation ( 1.1 1|) can be turned into the following one 



Py 
A„ 



eB 

Py X 

C 

x y + ~ X p* 



(3.1) 



Next let us notice what we did in ( |2.9| )- (|2.10| ) in order to go from the Hamiltonian (H) 
in Ai to the new Hamiltonian (7i) in Ai: we just replaced ip a with the superfield $ a . 

1 



3 This is the "analog" of the transformation formula (q\p) = 



exp[ipq/h] between momenta 



and position eigenstates but here the role of position and momenta is taken respectively by p and \ p . 
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Let us do the same for the minimal coupling in Ai space for the case of a constant 
magnetic field given by eq. ( |1 . 7|) . It means the following: 



Py 



-Bx 



i 



$Py B& J 

c 



Expanding flOD in 9,9 and using ( |2.19| ) and ( |2.2(JD we get 



p y — i99\ 



p y — i99\ 



eB 



Comparing the terms with equal number of 9 and 9 



x + i99\ Px 



becomes 



(3.2) 
(3.3) 

(3.4) 



Py 



A„ + 



eB_ 

c 
eB 



Py X 

C 



(3.5) 



A, 



which are exactly the substitution rules (3J.) for the minimal coupling for 7i. So we 
can say that the superfield formalism provides a compact way, eq. ( |3.3| ), to write the 
complicated minimal coupling (|3.1|) . 

Let us now check if this compact way of expressing things via superfields is an acci- 
dent of the case of a constant magnetic field or if it holds in general. The Hamiltonian 
H of a free particle in a generic magnetic field is obtained via the minimal coupling 
and is 



H 



2ml 



Px 



~A X 



Py 



-A 



The associated H B of ( |2.18|) is then: 



m 



H B = —[P X --A X )+^(py--Ay)+—(p 



m 



+ P 



A 2 



(3.6) 



m 



-A, 

c 



-A 9H - A 9H A 9H 

Vx dx Py dy Pz dz 



m 



edA 



K, + 



edA, 



+- 



m 
1 



+— A„ + - 



c dx " 
edA 
c dx 1 



c dy 



e dA x 
c dz 



" \ i_ \ _j_ e &A 
c dy c dz 



edA z edA z edA 
A z H — A„ T H — A D „ + - A p 



(p* - 


~A X 




c 








— Ay 




C 






fp 2 - 


-A z 

c 



(3.7) 



m v c dx Px c dy Py 
So we see that this last expression can be obtained from the Ji of the free particle 

H = —X x p x + —\p y + —\ z Pz (3.8) 
m m m 
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via the simple substitution 



e . 

Vi — 'Pi~ -An 

°e_/ dA n A ( 3 - 9 ) 
V 



This is the minimal coupling for TC in a generic magnetic field generalizing the one in 
a constant magnetic field given by eq. (|3.5|) . We have indicated with A q . the spatial 
components of the gauge field. We want now to see if ( ^.9[ ) can be derived from the 
superfield generalization of the standard MC, i.e.: 

Pi — Pi-- c A qi (q) (3.10) 

I I 
$Pi „ $w_f^ 4 ($9) (3.11) 

Let us first notice that, if we neglect the c, c variables, we have that the following 
relation holds: 

r - 8 A 



This can be checked easily using eq. (|2.12| ). If we now expand ( |3.11| ) in 9, 6 and we 
neglect the c, c, what we get is: 

Pi - i96\ qi Pi - i66\ qi - e -{A qi + i69A qi ) (3.13) 

Comparing the terms above with the same number of 6, 6 we obtain 

e , 

Pi — > Pi- ~A qi 

X A + e vfA (3 - 14) 



which are exactly the relations (|3.9|) . So this proves that (|3 . 1 1| ) is the most compact 
way to write the minimal coupling for 7i. 

The reader may ask which is the physical reason why \ q . should be changed, like we 
do in eq. ( |3.9p , when we turn on the magnetic field. To answer this question we have 
to analyze the issue of the gauge invariance of the system. Let us remember that the 
Lagrangian associated to the H of eq. ( |3.6|) was 

1 e 
L = -m(x 2 +f + z 2 ) + ~(xA x + jjA y + zA z ) (3.15) 
— c 
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where 



x 



— (Px 

m 

m c 

~{Pz - -A 
m c 



-a. 

c 

~A„ 



(3.16) 



The velocities which appear above are measurable quantities and so they must be 
gauge-invariant. Since A q . transform under a gauge transformation as A q . + d qi a(q) 
with a(q) an arbitrary function, the pi must transform as 



Pi — > Pi + -dg z a(q) 
c 



(3.17) 



The Hamiltonian H, being a combination of gauge-invariant quantities like Pi A qi , 

is gauge invariant while the L of eq. ( 3.15 ) changes by a total derivative. As the 
Hamiltonian H is basically the energy of the system it must be gauge invariant. Let 
us now ask ourselves how A should change under a gauge transformation. If we adopt 
the compact notation (|3.11| ) we used for the minimal coupling of 7i, it seems natural 
that the superfield analog of the gauge transformation ( |3.17| ) should be 



and 



$p* — > $ft + _[0 ft a]($*) 



A qi (&) — A v (&) + [d 9t a](& 



(3.18) 



(3.19) 



where by [d qi a]($> q ) we mean that we insert $ 9 in place of q in the function that we 
obtain by making the derivative of a(q) with respect to Expanding ( [3.18] ) in 6, 9 we 
get (|3.17|) as first component and the following one as last component 



x qi + 1 J dede - c [d qi a\m = K - - c KAAMq) 



Similarly expanding (|3.20|) in 6, 9 we get the usual transformation A q . 
as first component and the following one as last component: 



A n 



A' qi = A qi + d qi a(q, X p ) 



where 



a(q,X p ) 



(3.20) 

+ d qi a 

(3.21) 
(3.22) 



ic f — 

It is then easy to see that the combination X m / d8d8A q .(Q 9 ) = X q 

which is the X qi components of eq. ( |3.11| ), is gauge invariant if we gauge transform 
X qi as in ( ^20|) and A qi ($ q ) as in (gig ). So the RHS of Q) are gauge invariant 
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quantities and, as a consequence, also the Ti B of eq. ( |3.7| ) is gauge invariant because it 
is built out of the combinations (|3.9| ). Of course all this is very formal and stems from 
the extension of the standard gauge transformations to the superfields, ( |3.18| )-( |3TT9|) . 
While H, being the energy, must be gauge invariant there is apparently no reason why 
Ti should be gauge invariant. Similarly, while we know that pi should change under 
a gauge transformation like in ( |3.17 ) in order to make gauge invariant the velocities 



( |3.16| ) which are observables, there is apparently no physical reason why X qi should 
change under a gauge transformation as in (|3.20| ). Actually there is a physical reason 
and it is the following. As the velocities ( |3.16| ) are gauge invariant then their evolution 
has to be gauge invariant too. The evolution can occur via H and the standard Poisson 
brackets 

{<p a ,<p b } pb = u ab (3.23) 

or via Ti and some extended Poisson brackets { , } ep b which were introduced in ref. 
Q (for a brief review see appendix A). In this extended formalism the Hamiltonian Ti 
( |3.7| ) can be written, using the notation ( |3.12| ), in the following compact way: 



n- l -Y.{K- e -A){ P ,- e -A,) (3.24) 
The gauge invariant velocities 

v^ = -( Pi --A q ) (3.25) 

m\ c J 

evolve via the extended Poisson brackets in the following manner: 



v. 



{v gi ,Ti} epb (3.26) 



For example the result for v x is 



g 

Vx = — (B z v y - B y v z ) (3.27) 
mc 

If we now use a different gauge and just change the gauge field A qi and pi, but not X qi , 
we would get as new Ti: 

h'4e(a„-^)( p ;-H) < 3 ' 28 > 

where with A^,p\ and A' q . we indicate the quantities gauge transformed according to 

(Q)-ra- 1 



The evolution of the velocity v x , via the gauge transformed Ti', would turned out to 
be: 

v x = {v x ,H'}e P b = — (B z v y - B y v z ) H [(<%ai)v x + (d y d x a)v y + (d z d x a)v z ) (3.29) 
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So we notice that the evolution is not anymore gauge invariant because it depends on 
the gauge parameters a which appear on the RHS of Q3.29Q . This is absurd because 



the velocities are gauge invariant quantities and so their evolution must maintain their 
gauge invariance. This lack of gauge invariance is the price we would have paid by not 
allowing \ qi to change under a gauge transformation or, equivalently, by not allowing 



A to enter Ti via the MC combination of eq. (3J3). Note that the X Pi , differently than 
the X qi , do not have to be changed at all in order to maintain the gauge invariance of 
the system. This answers one question we raised in the introduction and concludes the 
analysis of the issue of the gauge invariance. 

Up to now we have regarded Ti as a function endowed with its own extended Poisson 
brackets. We want now to proceed to analyze the same issue of gauge invariance when 
we turn TC into an operator Ti, like we did in formula ( [2. 27\ ). Let us first briefly 
review what happens in QM, following ref. ||. In classical mechanics the gauge 
transformations leave q invariant but change pi as follows 

ti=Vi+-d qi <*(q) (3.30) 

and, as a consequence, the Poisson brackets are left invariant under these transforma- 
tions. 

tii,Pj} = UhPj} = % (3-31) 
Then the quantization rules for the transformed variables are 

tiuP'j} = $ij — ► [qi,pj] = ihdij (3.32) 

d 

This implies that p\ can be realized operatorially like the original i.e. p'- = —ih— — . 

J J dqj 

The quantum Hamiltonian is then 

H = (-^^- - / 2m ( 3 - 33 ) 

As the p[ have been realized as before, if we do a gauge transformation we get that 
only Ai change in H and the new Hamiltonian is: 



H' 



2 



2m (3.34) 



It is easy to check that one can pass from H to H' via a unitary transformation 

H 1 = UHU- 1 (3.35) 

where U = exp( i—a(q)) . So, differently than for the function H of CM, the Hamil- 

V en J 

tonian is not gauge invariant, but what is important is that the expectation values are 
gauge invariant. In fact, if the H transforms as (|3.35| ), the states change as 

W) = U\ij>) (3.36) 
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In the (q\ representation this becomes 



ip'(q) = exp[ *^«(?) mo) 



(3.37) 



which is the usual transformation by a phase under gauge transformations. We can 
notice that the expectation values of (ip'\p[\ip') and (ip\pi\ip) are related exactly as the 
classical momenta in ( |3.30| ): 



(3.38) 



Let us now turn to the KvN operatorial theory and check how the gauge transfor- 
mations are implemented. At the operatorial level we have to construct everything so 
that the expectation values of 



(ip\pi AM 1 ) and of 



(3.39) 



would be gauge invariant. Let us start by noticing that the commutation relations 
[(p a , Aft] = i5% are the operatorial counterpart of the extended Poisson brackets 
{(p a , \b\epb = 8% an d ^he S au S e transformed coordinates ip' a , X' b under fl3.17| ) and fl3.2(J| ) 
have the same epb {y? /a , X' b } ep b = 8% as the original variables. So we expect that also the 
associated commutators among the gauge transformed operators would be the same as 
the original one: 

W a X\=< (3-40) 

This means that we can represent the (p' a and X' b in the same manner as the (p a and 
Aft. As a consequence the gauge transformed version of the quantities in ( |3.39| ) is: 



(m--/ qi -- c [d qi a)(q)m 

m\ qi -lA qi -~[d q MQ,KW) 

Note that, via the introduction of the following operator 



(3.41) 



U = expl - i-X Pi [d qi a](q) 



we can write the following transformations: 

Pi ~ -An - -[d qt a](q) = U 



Pi 



-A 



u- 1 



A —-A 



■[d qi a](q,X p ) = U 



X —-A 



u- 1 



(3.42) 



(3.43) 



This implies that ( |3.41 ) will be gauge invariant provided we transform the states as 
follows 



W) = U\i)) = exp\ -i-X Pi [d qi a](q) 



(3.44) 
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Let us now represent this transformation law on the two basis given by ( p. 31 ) and 
In the basis (|2~3^ ) we have from flOp 



exp^-i-X Pi d qi a(q) \(qX 
exp^i-a^jip(q, X p ) 



(3.45) 



So in this basis the gauge transformation is just the multiplication by a local phase 
factor a in the space (q, X p ) where a has been defined in ( |3.22| ). 

Now let us represent ( |3.44|) in the {q,p\ basis (|2.31|) and let us make use of the 
transformation formula (|2.34|). What we get is: 



4>'(q,p) = {QPW)= dq'd\'{qp\q'\'){q'\'\U\ij) 



dXL 



=exp[ipX' p ]exp 



i-a(q,X') 
. c 



(3.46) 



Inserting a further completeness we obtain 



ip'(q,p) 



dX p dq'dp'exp[ipX' p }exp i-a(q,X' p ) (qX' p \q'p')(q'p 



dX p dp' 
2tt 



exp 



iX' p (p — p da) 



J dp'5^p — p' — -dc^Jip(q,p') 



ip(q,p') 
ip(q,p 



-da 

c 



(3.47) 



So in the (q,p) representation of our Hilbert space the gauge transformations are not 
implemented by the multiplication by a local phase, like in the (q, X p ) representation, 
but by just a shift in the argument p of the wave function. It is easy to show, as we 
will do in the appendix B, that the Liouville eq. ( |1.2|) is invariant in form under the 



gauge transformations. In appendix C we will show that the phase expi i-a ) of ( p. 45 ) 



is exactly the one that can "pass through''^ the Hamiltonian 7i of KvN if we change 



the gauge field as we do in ( 3.19 ). In that appendix we shall also explore which gauge 
fields have to be inserted in H to allow for a general^] phase exp[ia(q, X p )} to "pass 
through" H. 



6 By "pass through" we mean a procedure explained in appendix C. 
7 By "general" we mean one not of the form expi i-a 
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4 Landau Problem 



In this section we will make a first application of the minimal coupling scheme that we 
have developed previously for the KvN formalism. This first application is the Landau 
problem. We will first review it in quantum mechanics and then turn to the classical 
operatorial version of KvN. The Landau problem is concerned with the dynamics of 
a particle under a constant magnetic field directed along z. We make the following 
choice for the gauge potential: 



A x — 0, Ay — Bx, A z 







The Schrodinger Hamiltonian is then 



H 



1 

2m 



Px + [Py- —X 



Pz 



(4.1) 



(4.2) 



As p y ,p z commute with H we can diagonalize all these three operators simultaneously. 
The eigenfunctions will then be labelled by the eigenvalues of H, i.e. E, and by those 
of p y ,p z which are Py and p° z . Their form will be 



2nh 



exp 



Ti 



{ply + p\z) 



The stationary eigenvalue problem 

Hip BtP o >p o = E?p E j> iP o 
leads to the following differential equation for ijj(x) 



if>(x) (4.3) 



(4.4) 



2m^ ^ ^ 2m 
Indicating with E t the quantity 



6 Bx \ ^ 

— — J ip(x) 



p 



.o 2 



Et = E — 

2m 



and making a change of variables from x to x', with x' = p G y 
into the following one 



eBx 



x) (4.5) 



(4.6) 



, eq. ( fL5| ) is turned 



ti 2 I 
_ —f'( x ') + — 

2m 2m 



1 2 



eB 



x' 2 ip(x') 



eB 



(4.7) 



We can immediately notice that this is like an harmonic oscillator eigenvalue problem 
with the frequency replaced by 

lj = (4.8) 
eBm v ; 
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and with the energy replaced by 
harmonic oscillator problem: 



eB 



E t . This quantity is discretized like in the 



1 2 



eB 



(4.9) 



Combining 



and ( |4.8|) we get from ( |4.9|) 
ehB ( 



mc 



\ V 2m 



(4.10) 



So the eigenf unctions ( |4.3| ) can be labelled by the quantum numbers (n,p°,p°): i.e. 
V'n.pQ.pO- Note that these wave functions are degenerate because all those with different 
values of p° y have the same E n ^ p o . So there is an infinite degeneracy. 

Let us now analyze the same problem at the classical level using the operatorial 
formalism of KvN. Using ( |3.24j ) and the gauge choice ( |4.1| ), the 7i is 



n 



-Kpx 



A,, - -Ay) [Py - -A, 



H — Kp z 



(4.11) 



We will now turn 7i into an operator TC using the "mixed" representation ( [2.30|) and 
what we get is: 



n 



1 d d 

m dx dX r . 



1 / . d eB , 
m\oy c 



d eB \ 1 8 d , A , n . 

7Tx x +-^-^— ( 412 ) 

OAr>„ c J moz oXr 



Let us now diagonalize this operator. The reason to do that is because in the KvN 



theory the equation to solve is (|1.2|) which, because of (|2.27|) , can be written as 

id t i\) = ity (4.13) 
So, like for the Schrodinger equation, one should first diagonalize TC 

Hi>~ E = E^~ (4.14) 
and then write a generic wave function as 



-iEt 



(4.15) 



where the C7g are derived from the expansion of the initial ip on the ife. We want 
to underline that the E that appear in eqs. ( |4.14| )-( [4.15| ) have nothing to do with the 
physical energy of the system. They are simply the possible eigenvalues of the evolution 
operator H and, using them and the associated eigenfunctions, we can deduce the 
evolution of the ip like it is done in formula ( J4.15 ) . For more details about this and on 
the manner to reconstruct the standard deterministic motion of CM see ref. 01. 
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Let us now turn to ([4.12 ) and diagonalize it like in ( |4.14| ). Note that the operators: 

d d d d 

— and i— — commute with Ti and so we can diagonalize simulta- 



dy dX Pv dz 



OXr. 



neously these five operators. The generic eigenfunction (in the mixed representation 
q, Xp) has the form 



A ^ = j2^y exp \- iX °y y ~ iX PyP°y\ ex Pi iX °z z ~ iX PzP°z\Tp( x ' V) ( 4 - 16 ) 



d . d 



d . d 



re- 



where A„, A" and pi are eigenvalues of the operators —i^r, - , ■ 

y v dy dX Py dz dX Pz 

spectively. We see the similarity with the quantum case except for the fact that the 

dimension of the space is double. Inserting ( [4.16|) in (|4.14j ) we get the equation 



1 d d 



eB 



-—— + -[Xl + —X 

m ox X Px m 



Via the new quantity 



we can rewrite eq. ( |4.17|) as 



7Ar 



eB 



Py ~ ~~ x ) + -ASP? 



in 



e + = e- -xy z 

m 



ij(x,X Px )=Eij(x,X Px ) (4.17) 



(4.18) 



1 d d 
m dx dX n 



ni 



eB 



+ -[K + —K*[Pv-—x 



eB 



i;(x,X Px )=E + i;(x,X p J (4.19) 



Doing now the following change of variables 



/ c o 

X = X ~ ^B Py 

K* = V + ~^ X y 



we can rewrite eq. ( |4.19D as 



l__d d_ ]_feB\ 2 x , i 

m dx' dXL m \ c J 



•Pm" 



^(x',X' Px ) = E + ^(x',X' Px ) 



(4.20) 



(4.21) 



The dimensions of the various quantities are such that we can write the above equation 
as 

' ' <} (> '■' ' ' - - ' ' ><< - (4.22) 



m dx' dX 



muo A„ x 



P-r 



ij(x',X' p J = E + ^(x',X' Px ) 



eB 



where uo = — has the dimension of an angular velocity and it is related to the well- 

mc 

known Larmor frequency of rotation of a particle in a constant magnetic field. Eq. 
([4.22|) is the KvN eigenvalue equation for an harmonic oscillator which is studied in 
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details in appendix D which we advice the reader to go through before going on. The 
spectrum of the 7i for the harmonic oscillator (see appendix D) is given byQ 

E° + s ° = Noj } N — ■ ■ ■ , —2, —1, 0, 1, 2, • • • (4.23) 



So, using ( f4.18| ), the final spectrum of eq. ( |4.17|) is 

E = N(—) + -\y z (4.24) 
\mcj m 

and the wave functions are any linear combination of the following eigenfunctions (see 
appendix D): 

^7v,n,AO, p o iA o iP o = ^-exp[iX° y y - iX Pu p y }exp[iX° z z - i\ Pz p° z ] ■ < sc (Z+)?/&(Z_) (4.25) 

where ip^ sc are the eigenfunctions of the quantum 1-dim harmonic oscillator^ with h 
replaced by an arbitrary quantity A which has the dimension of an action and with 
Z + , Z_ defined as : 

, x' + AX' x'-AX' 
We see from the form of the wave functions in ( |4.25| ) that the degeneracy is much more 



than in the quantum case. Not only the eigenfunctions with different values of p y have 
the same E, but the same happens for those eigenfunctions with different values^ of 
Xy, n and of X° z ,p° z . So it is a much more wider degeneracy than in the quantum case. 
The reader may wonder why this happens. We feel that this may be due to the fact 
that the "wave functions" in the KvN formalism have a number of variables {q,X p ) 
that is double than in QM. 

The reader may be puzzled by the presence in the eigenfunctions ( (4.25 ) of the arbi- 



trary quantity A which did not appear in the original 7i. This is not a problem. In 
fact once we are given an initial wave function ip(q, X p ) not depending on A, we will 



expand it on the basis ft4.25|) as: 



^(<L \>)= 53 (^n.Ag^A^g) (V , AT,n,A>°,A°,p°) (4.27) 

where Cjv, n - are coefficients which will depend on A and this dependence will com- 
pensate the one contained in the ipjsr^-- The evolution in t will not reintroduce the 
dependence on A because the eigenvalues of H do not depend on A as shown in ap- 
pendix D. Solving a classical system via its KvN states is like working in the Schrodinger 
picture. We could also work out the analog of the Heisenberg picture and this is done 
for the Landau problem in appendix E. 

8 Note that we have a discretization phenomenon even at the classical level for the eigenvalues of 
the Liouvillian. This is related to the requirement of single valuedness of the KvN states as explained 
in appendix D. 

9 n = -N, -N + 1, —N + 2, ■ ■ ■ if N is negative and n = 0, 1, 2, • • • if N is positive. 
10 with the only constraint that the product X z p z be the same. 
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5 Aharonov-Bohm Phenomenon 



The second application of the MC that we will study is the well-known Aharonov-Bohm 
(AB) effect, 0. This is a phenomenon which proves that the QM wave functions are 
changed by the presence of a gauge potential even in regions where the magnetic 
field associated to this gauge potential is zero. The change in the wave functions can 
be detected by an interference experiment. The classical motion instead feels only 
the magnetic field and not the gauge potential. In this section we will study this 
phenomenon not by looking at wave functions but at the spectrum of respectively the 
Schrodinger operator H and the classical KvN Liouville operator 7i. The geometrical 
set up that we will use for the AB effect is illustrated in Figure [l| and it has been 
suggested in the book of Sakurai 0. Basically we have two infinitely long cylinders 
one inside the other. We will study the Schrodinger operator H in the region in 
between the two cylinders. We will show, as Sakurai H indicated, that the spectrum of 
H changes once we turn on the magnetic field inside the smaller cylinder. In the region 
in between the two cylinders the magnetic field is zero because the smaller cylinder 
shields completely the magnetic field. So we are exactly in an AB configuration: zero 
magnetic field and non-zero gauge potential. Using the same geometrical configuration, 
we will study the spectrum of the KvN-Liouville operator 7i and we will prove that it 
does not change once we turn on the magnetic field differently than what happens in 
QM. We feel that this, in the framework of the operatorial formulation of CM, is the 
best mathematical proof that there is no AB effect in CM. 

Let us now study the Schrodinger operator in the geometrical set up of Figure [1] and 
let us do it at first without magnetic field. The Schrodinger operator in cylindrical 
coordinates 

x — pcosip 

y = p simp (5.1) 
z = z 



is for a free particle: 



H = -—(— + -— + -— + — ) (V?) 
2p J \dp 2 pdp p 2 dip 2 dz 2 } 



and the eigenvalues equation is 

h 2 ( d 2 Id Id 2 d 2 \ 

-tAw + Wp + 7W + d^r ip ' z) = E ^ p ' ^ z) (5 - 3) 

d d — 
As the operators — — and — commute with H we can diagonalize these three operators 

dip dz 

simultaneously and have 



1 \ip\z 



^(p, <P, z) = ^exp 



exp[irrup]R{p) (5.4) 
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where p Q z is a fixed value and m is an integer. Inserting (5A) in ( |5.3|) we get the following 
equation for R(p) 

R"(p) + ^+(-s-^r)R(p)=0 (5.5) 



p \ p 
where 

- _ 2fiE p° z 



(5.6) 



h 2 h 2 

If s = 1 then eq. ( |5.5|) would be the well-known Bessel equation 0. We can get it by 
using the new variables r = v^P- I n " r " ec i- (|5.5|) is 

d 2 R ldR / m 2 \ n 

+ -^-+ I--5- fl = (5.7) 



As we want the wave function to be confined between the two cylinders, the boundary 
conditions should be 

^(^a) = R{sfsb) = (5.8) 

where a and 6 are the radii of respectively the smaller and the larger cylinder, see 
Figure [I]. The general solution of eq. ( |5.7| ), with m integer, is given by the linear 
combination of the Bessel functions of the first and second kind which are §J : 



n=0 n ' +Tn+i) 

Y m (r) = lim-[J m+£ (r) - (-l) m J_ m _ e (r)] (5.9) 

The general solution will then be 

R(V^p) = AJ m {^p) + BY m {y^p) (5.10) 



Imposing the boundary conditions (|5.8j) we will get the spectrum of the system. In 
order to simplify things we will consider the limiting case in which the radius of the 
internal cylinder a goes to zero. In this case the boundary conditions ( |5.8|) become 

R(0)=0, R(V^b) = (5.11) 

It is well known H that the Bessel functions of the second kind Y are singular in the 
origin, so we restrict ourselves to solutions ( |5.10| ) of the form: 

R(V§p) = AJ m (V^p) (5.12) 

With this choice the first of the boundary conditions ( |5.11| ) is automatically satisfied 
because J m (0) = for m > 1. We have to satisfy only the second one of the 
conditions ( |5.11| ) which implies: 

J m (Vsb) = (5.13) 
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This relation tells us that we have to look for the zeros of the above Bessel functions. 
Let us call them a^^ m where m indicates to which Bessel function we refer to and k 
labels the various zeros of the m-Bessel function in increasing order, so k = 1, 2, 3, • • •. 
The solutions of eq. fl5.13|) can then be formally written as 

\^b = a k)m (5.14) 

Replacing s in the equation above with its expression ( |5.6|) , we get that 



E, 



(5.15) 



These are the energy levels. If we choose m — 1 and the second zero (k = 2) which is 
«2,i ~ 3.83, we get 



(5.16) 



fib 2 2fi 
see Figure |]. 

Let us now turn on the magnetic field |5| which we want to be zero everywhere except 
for the B z component on the line x 2 + y 2 = and with a fixed flux $ s . A choice of 
the gauge potential || is 



A, 



A, 



2n(x 2 + y 2 ) 

X<& B 



y 2ir(x 2 + y 2 ) 
. A z = 

Turning to cylindrical coordinates we have that ( J5.17|) is equivalent to 

A p = 0, Ap = — — , A z = 
2np 

So with the above choice of gauge the minimal coupling affects only the — ih 



(5.17) 



(5.1i 



d_ 

dip 



— ih 



d 



-ih 



d e$ ; 



dcp dip c 27r 

The associated Schrodinger operator (|5.2|) is then: 



(5.19) 



H 



2/i 



h 2 r d 2 id i 



dp 2 pdp p 2 \ dip ch 



d 



/e 



+ 



d_ 

d? 



2 i 



(5.20) 



and the eigenvalue equation is 



-h 2 



2/i 



d 



+ 



Id Id 2 
+ 



lie 



1 d 



dp 2 pdp p 2 dip 2 ch B p 2 dip c 2 h 2 p 2 ' dz 2 



e 2 $| d 2 
+ 



i/> = Ei/>(p,<p,z) (5.21) 
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Like for 



we can choose solutions of the form 



^{p, ^> z ) = ^-exp 



ip° z z 



exp[irmp]R(p) 



(5.22) 



which, inserted in (|5.21|) , give the following differential equation for R(p) 

2 jS 2 ■ 

1 I m. 

2/i 



h 2 \R" 1 R' 1 
-R + JR-7 [V> 



ch 



p\_ 

h 2 



E 



(5.23) 



e $ s 2fiE p 

Indicating with a = — — and with s = — n V, eq. (15.231) can be written as 

ch h h — v 



P 



[m — a) 



s — 



R{p) = 



(5.24) 



If we compare the previous equation with ( |5.5| ) we notice that having turned on the 
magnetic field has only shifted m — > m B = m — a. Here m B will not be anymore an 
integer but a real number. Doing the same change of variables as before, r = \^fp, we 
can transform ( |5.24| ) into 



d 2 R 1 OR 



m. 



-8^ + rfr + { 1 -^) R = 



(5.25) 



For this equation with m B real there are two linearly independent solutions which are 
two Bessel functions of the first kind with opposite indeces: 



Jm B {r) 



J-m B (r) 



!)»(§ 



r\ m B 00 

y — 

^ n! r(n + m 



,2h 



1) 



2n 



^ n!r(-n + m s + 1) 



(5.26) 



One immediately notices that J ms (0) = while J_ ms (0) diverges. As before we must 
have as boundary conditions J ms (0) = and so the general solution of (|5.25|) is 



R{r) 



AJ mB (t) 



The other boundary condition gives the following relation 

Jm B {V^b) =0 

from which, as before, we can derive the energy levels 

zp _ f- a k,m B v\_ 
k ' mB ~ 2pb* + 2/i 
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(5.27) 



(5.28) 



(5.29) 



If $ B is such to give, for example, a = 0.1, then we will have to consider the Bessel 
functions J m _o.i- The second zero of J .g, analog to the second one of J\ that we 
considered before, is a 2 ,o.9 — 3.70, see Figure 0. Inserting this value in (|5.29|) we get 



2,0.9 



2,1 



(5.30) 



The last inequality indicates that £"2,0.9 is smaller than the corresponding level £2,1 of 
the case without magnetic field calculated in ( |5.16 ). So this is a clear indication, like it 
was suggested in ref. ||, that the presence of a gauge potential modifies the spectrum 
of the Schrodinger operator even if the wave function is restricted to an area with zero 
magnetic field. 

Let us now perform the same analysis in the classical case using the KvN operator 

7i. We have to be careful here because, we have to go to cylindrical coordinates both 

d d 

for the q and the p variables and for their derivatives — , — which enter Tt. In the 

oq op 

case without magnetic field the Lagrangian is 



11 1 
L = -pp 2 + -ppV + -pz 2 



(5.31) 



and so the momenta conjugate to p, <p, z are: 

V P = PP 

Pv = ppV 

Pz = H'z 

The relations between p x ,P y ,Pz and Pp-.Plp-.Pz can be easily worked: 

P 9 . 

p x = hx = fipcosif — Hpipsvnip = p p cos(p -simp 

P 
v 

p v = uy = upsirup + ppipcosu) = PpSirup H — -cosip 

P 



(5.32) 



(5.33) 



I Pz 



Pz 



As the z and p z are the same in the two coordinate systems we can summarize the 
basic transformations in the following set 



x = pcostp 
y = psirup 

Px = PpCOSLf 



2V 
P 

Pv 



-sirup 



(5.34) 



Pv = PpSirup H — -cosip 
P 



We should note that the transformations of the momenta p x ,Py are not just functions of 
the new momenta p p ,p v but also of the cylindrical coordinates <p, p. Let us remember 
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d d 

that 7i in cartesian coordinates contains the derivatives — , — and so we should check 

dq dp 

how they are related to the derivatives in cylindrical coordinates. Using (|5.34 ) it is a 
long but easy calculation to show that: 



sirup cosip/p p^cosip/p 2 p v /p ■ sirup — p p cos<p 
cosp —psinip 

sirup pcosip 



( d/dp \ 
d/dp 
d/dp p 
V d/dp v ) 
(5.35) 

Equipped with these transformations we can then easily transform 7i from cartesian 
coordinates to cylindrical ones. In the case of a free particle we get 



/ d/dx \ 


( 


d/dy 




d/dp x 




V d/dpy ) 


\ 



rT, -Px d .p y d .p z d 

p dx p dy p dz 

d i pm d i d 

Pz 



iP% d 



pF p dp p p 2 dp p^ z dz p p 3 dp p 



(5.36) 



where in the second step we have used ( |5.34| ) and ( |5.35| ). Next we will turn on the 
magnetic field whose gauge potential is in (|5.17| ). The minimal coupling for 7i is given 
in ( |3.24j ) and we need it to build the A qi . For our potential these A qi are 



A^ 



$ B xy $ a x 2 - y 2 

Px n (x 2 + y 2 ) 2 Pv 2vr (x 2 + y 2 ) 2 

$ s x 2 - y 2 $ fl xy 

An., — ; — Ti ^tt: + Ar 



(5.37) 



y Pv 2ir (x 2 + y 2 ) 2 
A z = 



7T (x 2 + y 



2\2 



In the expression above we have now to turn the \ Px , X Py into operators, like in Q2.26| ), 
and next we have to change everything into cylindrical coordinates using (|5.34j) and 
( p.35| ). The result is 



A x 

A = - 
A z = 



1 2tt p 2 
1 2tt p 2 



strip 







pcosp- 







dp p dp^ 

d d 

cosp- h psmp- — 

op p op v 



(5.31 



Inserting (|5.37|) - fl5T38D into ( |3.24|) and turning the X q . into operators, after a long but 
trivial calculation, we get 



T-Ca = PpT, n[Pip 

p Op pp l ( 



e$ R \ d 



i d 

Pz 



2irc J dp p dz pp 



P v 



e$, 



d 



2ixc J dpp 



(5.39) 
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Notice that we could obtain this Ti A from the Ti of eq. ( |5.36|) by just doing the 
replacement 

p * - ** - Wc (5 - 40) 

The eigenvalue equation in the free case 

Ti ip(p, p p , <p, p^, z, p z ) = E ip(p, p p , <p, p v , z, p z ) (5.41) 

d d 

can be solved by noticing that Ti commutes with the four operators — i-^— ^tt- iVviPz- 

dip oz 

So these five operators can be diagonalized simultaneously and the solution will have 
the form 

ip = —R(p,pp)5(p v - p°)5(p z - p z )exp(imp)exp(i\° z z) (5.42) 

d 

where n is an integer and p^,,p° z , X° z are the eigenvalues of p v ,p z , —i—. Inserting ( |5.42|) 



in ( |5.41| ) we get the following equation for R(p, p p ) 

-PpJT + 3^1 7^ + — -E)R(p, Pp ) = 5.43 

p dp pp 2 pp A v dpp p J 



We showed before that we could turn Ti into the Ti A of ( p.39| ) by just doing the 
substitution ( 5.40|) p v —> p 9 — — — . It is then clear that we can turn also the solutions 
( |5.42j ) into the solutions of the eigenvalue equation 

H A ip A = E a i/j a (5.44) 

by just doing the substitution (|5.40|) into ( |5.42|) . The result is: 

ip A = -^Ra{p,P p )s(p v -pI - ( ^)6(p z - p° z )exp{inip)exp(i\° z z) (5.45) 

Once we insert this into eq. ( |5.44|) we will get for R A (p,p p ) the following equation: 

.«» + ^ - J-rf JL + M _ eY RaM = (5.46) 
p op pp z pp A r opp p J 

So eq. ( |5.46j ) is the same as the free one ( |5.43| ) and as a consequence the spectrum E A 
is the same as the E of (|5.43| ). This is the proof that the spectrum of the Liouville 
operator is not changed by the presence of the gauge potential. Of course the two 
eigenfunctions which have the same eigenvalues E = E A are different because they are 
labelled by different eigenvalues of the operator p v . In fact the eigenfunction ift of eq. 
( |5.42j ) has eigenvalue p° for the operator p^ while the eigenfunction ip A of eq. Q5.45 ) 

has eigenvalue p° H . So the two eigenfunctions are related by a shift in one of their 

^ 2-7TC 
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"classical" numbers p Q . The difference with the quantum case is that the corresponding 
eqs. ( |5.5|) and Q5.24Q cannot be turned one into the other like we did in the KvN case. 
This is so because in (|5.5| ) m is an integer and not a continuous real eigenvalue like 
p° is in the KvN case. Also in classical mechanics we had an integer eigenvalue, n for 
d 

but we managed to down- load on p° and not on n, the difference between the 

dtp ' v 

free and the interacting case. 

The reader may object that, even if the classical spectrum is the same in the two 
cases, the eigenfunctions are different and then the evolution may lead to different 
results. Actually it is not so because, as we see from Q4.15Q , we have to integrate over 
all the possible eigenvalues which label the eigenfunctions. In our case the different 
eigenfunctions (|5.42|) - (|5.45|) have only the "classical" number p° shifted. Since p° can 
assume every real number, when in ( |4.15| ) we integrate over all the p° a shift in them 
has no effect on the final result. 

We feel that this proof that the spectrum of the classical KvN operator is unchanged 
by the presence of the gauge potential, while the spectrum of the Schrodinger operator 
is changed, is the most convincing proof of the AB phenomenon. 



6 Conclusions 

In this paper we have studied which is the minimal coupling procedure for the KvN op- 
eratorial approach to CM. We have shown that the MC involves not only the momenta 
but also their derivatives. We managed to encapsulate these two MC into a single one 
using the concept of superfield. We have then applied this technique to the Landau 
problem and to the Aharonov-Bohm phenomenon. In the first case (Landau problem) 
we showed that in the KvN formalism there is a sort of discretization phenomenon in 
the eigenvalues of the evolution operator. Moreover we proved that there are many 
more degeneracies in the classical than in the quantum case. For the second problem 
(the Aharonov-Bohm one) we showed that at the quantum level there is a change in 
the spectrum of the Schrodinger operator once the gauge potential is present while 
there is no change in the spectrum of the classical KvN operator. We feel this is the 
most convincing proof of the AB effect. The paper contains also a complete analysis 
of the issue of gauge invariance in the Hilbert space of KvN. 

Having now all the tools to write down, in the KvN formalism, the interaction 
between a particle and a gauge field, what we should do next is to see how the gauge 
fields interact among themselves in the KvN operatorial approach. This has already 
been started and work is in progress on it ||10| . 
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Appendices 



A Appendix 

In this appendix we will briefly review the extended canonical formalism associated to 
the Hamiltonian 7i of eq. ( |2.8| ). From the Lagrangian ( |2.7| ) one could derive the equa- 
tions of motion for the 8n-variables ((p a , X a , c a , c a ) by the simple variational principle. 
These equations are: 

<p a -cu ab d b H = (A.l) 
[5 a b d t - uj ac d c d b H}c b = (A.2) 
5 a b d t c a + c a oo ac d c d b H = (A.3) 
[5 a b d t + u ac d c d b H]\ a = -ic a u ac d c d d d b Hc d (A.4) 

We could ask ourselves if these same equations could be derived from the Hamiltonian 
7i. The answer is yes. If we introduce the following extended Poisson brackets structure 
in the extended space (ip a , X a , c a , c a ) 

W a Ab} epb = S a b (A.5) 
{c b ,c a } epb = -i5 a b (A.6) 

(while all the other brackets are zero) we get that the equations of motion (|A.1|) - (|A.4|) 
can be derived as 

^ = {0,H} epb (A.7) 

where O is any of the variables (ip a , X a ,c a , c a ) or any function of them. More details 
can be found in ref. H. 



B Appendix 

In this appendix we will prove that the Liouville eq. ( |1.2| ) is invariant under the gauge 



transformations. Let us write (|1.2j ) in the abstract form 



ij t \m)=n\m) (b.i) 

Let us now do a gauge transformation by a parameter a. The new ket will be 

W{t)) = U\m) (B.2) 

where U is the expression in formula ( |3.42| ). We will prove that this state satisfies the 
following equation 

*j t W(t)) = <H>W{t)) (b.3) 
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where TC' is the operator obtained from 7i by doing a gauge transformation 



H' = -[\ n 

m 



—X n ,„d n ,„ A„ 
c 



Pi 



-A' 



e\ v .d q .® 



(B.4) 



where the A' is the gauge transformed vector potential and $' is the gauge trans- 
formed scalar potential. ( |B.4p is the 7i in the gauge (A', $') and it is derived from 



H 



2m 



+ e$'. Let us now evaluate the LHS of (P-3| ) 



ij t m)\m)] 

e da 



d 
dt 



U(t) 



e da 



+U(t)H\m) = -- orl^(*)> + U{t)HU- x W{t)) = —— + H\W{t)) (B.5) 



c dt 



c dt 



where 



n = unu 



-i 



(B.6) 



In the first steps above we have used the expression (|3.42| ) for U and formula ( p.22| ). 
The explicit expression for 7i is 



m 



« = - K + zKAAM [Pi - -AM - <AMq' 



(B.7) 



where 



g 

Ki = U K U ~ = K + -KAA^q) 



Pi = UpiU 1 =pi- -d qi a(q) 

q[ = UqiU~ l = qi 
X' = V\JJ~ l = \ n 



(B.8) 



Remembering that 



A' = A a . + d a a 

qi h% Hi 

$' = $ — 

c dt 



(B.9) 



we can rewrite the Ti. of ( B.7| ) as 



H 



— [x qi + -XpAAMq) + -\Ak A iM) )\Pi- ~ d qMq) - -AqM 



in 



-eA„ A $ 



m 



K + -\A k A ' )\Pi- -K ) - eX P Aj 



1 da 
c~dt 



e d „ 
cdt 



Using this result 



becomes 



ij t W)=H'W) 



;b.io) 



(B.ll) 



which is what we wanted to prove. 
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C Appendix 



We know that one of the effects of the introduction of the MC is that local phases 
multiplying the states can be absorbed by a gauge transformation of the gauge field. 
What we mean is the following: if 



ip'(q) = exp 



and 



then 



H 



1 

2m 



-ih 



d_ 

dq 



-Aq) 



H'ifj' = exp 



l di a{q) 



ip(q) 

-in— 
aq 

Hip 



-Aq) 



(C.l) 
(C.2) 



where H' is obtained from H by replacing A with its gauge transformed A' = A + 



(C.3) 
da 



dq 



Now if we do the minimal coupling at the level of 7i, like in ( |3.24| ), the phase that can 

e 
i-a 



"pass through", like in ( |C.3| ), is exp 



c 



exp 



with a given in ( |3.22| ): 

Hip 



.e _ 
i-a 
c 



(C.4) 



The proof goes as follows. Let us use the mixed representation ( |2.30|) for the H asso- 
ciated to (IH): 



1 

m 



d e, dA\f d 



Av) 



and let us then transform TH, into a ii! where A — ► A' + 



da 
dq 



(C.5) 



In the same mixed 



representation ( |2.30D the wave function will be of the form ip(q, X p ) and we can also 
construct the following new state: 



ip'{q, A p ) = exp 



-i-X p d q a(q) 



exp 



.e _ 
i-a 

c 



(C.6) 



Equipped with these tools it is then a long but easy calculation to prove ( |C.4j ). To 
conclude we can say that the local phase transformations of the form ( |C.6| ) on the KvN 
states ip(q, X p ) are the classical counterpart of the local phase transformations ( C.l ) on 
the quantum Hilbert states ip(q). 

At this point a question which arises naturally is the following: if instead of the very 
particular phase transformation ( |C6[ ) on the KvN states we perform a general local 
phase transformation of the form 



4>'(q, X p ) = exp[ia(q, X p )]ijj(q, X p 



(C.7) 
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which gauge fields do we have to introduce in the 7i in order to absorb the phase like 
in ( C.4 )? The answer is the following: let us start from the 7i of the free particle 

1 d 2 



n 



m dqdX p 

and do a general phase transformation 

i>'(Q, K) = exp[ia(q, X p )]^(q, X p ) 
If we now perform in (p[H) the following MC 



— i- 



d_ 

dq 



■ 9 A 

- ? 7T + A 
oq 



._d_ 

'dxZ 



._d_ 



Ay 



(C.i 



(C.9) 



(CIO) 



where A q and A\ p are two gauge fields which transform as follows 



A' = A n 



da(q, \ 
dq 



A[ 



A\„ — 



da{q, \ 
dXr, 



then the new 7i, 



d 

n A = —[ -i— + A 



dq 



_d_ 



-Ay 



would satisfy the following relation 

H' A ip' = exp[ia(q, \ p )}H A ip 



(C.11) 



(C.12) 



(C.13) 



where H' A is the gauge transformed of H A via eq. ( |C.11| ). 

We should notice that ( |C.9| ) is a more general gauge transformation than the one in 
( p.6| ). To implement ( |C.9| ) we need two gauge fields A q , A Xp while for the ( |C.6| ) we 
could build everything from one field A(q), see ( C.5 ). The transformation ( C.6 ) is a 
particular case of ( |(J.9| ) and the same is true for the gauge fields which enter ( |C.5| ). In 
fact comparing ( |C.5| ) with ( |C . 1 2|) we see that the two general gauge fields of (|C.12|) 
should be of the following particular form in order to reproduce ( |C.5D : 



A n 



e OA 
c oq 



Ay 



-Mq) 



(C.14) 



It is interesting to notice that the interaction introduced by the two general gauge fields 
A q , A\ p appearing in (|C.12j ) is something new and not related to any magnetic field in 
g-space unless the fields have the particular form ( C7H ). We hope to come back in the 
future to the study of the forces generated by the general gauge potential of ( |C.12|) . 
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D Appendix 



In this appendix we want to build and diagonalize the Ti associated to an harmonic 
oscillator in one dimension. The Hamiltonian is 

H = — p 2 + -muj V (D.l) 
2ra 2 H y ' 

So the n (fngp is : 

H = \— - muj 2 \ p q (D.2) 
m 

which, in its operatorial version, using the representation fl2.26|), is: 



H = -i(^-mJ t q^-) (D.3) 
\mdq dp J 

This operator is very similar to the components of the angular momentum in standard 
quantum mechanics, so we can diagonalize it with similar techniques. First of all we 
can turn the standard phase space coordinates (q,p) into the following new ones (r, 6) 
defined as: 

p 

\pmujq = rcosO, —= = rsinO (D.4) 

The Hamiltonian becomes: 

H ■ r ■ gf^ r ® _j_ ^ ^ \ _)_ -y/ — o(^ T ® ® \ ■ ® (D 5) 

\fm \dq dr dq 39 / \dvdr dvd6/ 86 



V dq dr dqdO J \ dp dr dp d9 , 

The eigenfunctions of the previous operator are: 

xj){r, 9) = F{r)e~ im (D.6) 

and the associated eigenvalues are E = Nu. Imposing the single valuedness of the 
wave functions we get that 

il)(r, 6 + 2tt) = ^(r, 6) e~ 2ivN = 1 (D.7) 

which implies: G {0, ±1, ±2, . . .}. So the discretization of the eigenvalues of the 
Liouvillian is a direct consequence of the requirement of single valuedness of the KvN 
wave functions, requirement that was already present in Koopman's original paper Q. 

The proof of the discretization of the spectrum of the Liouville operator for an 
harmonic oscillator can be worked out also in the mixed representation ( [2.30| ) and 
( |2.33| ). In this representation the TC of eq. (|2.18|) becomes: 

# = -rnufiXtf (D.8) 

m dq dX p 
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Let us introduce the following new variables 



q + AA ? 
V2 



Z_ 



-A\ 
V2 



(D.9) 



where A is a constant which has the dimension of an action. In terms of these new 
variables the TC of (D.8I) can be written as 



H = -r 



A 2 d 2 



+ 



moo 



2mdZ 2 _ ' 2 
d 



1 

A 



A 2 d 2 moo 2 2 
2m~dZl + ^2~ + 
9 



(D.10) 



As indicated in the second step above, we notice that H is the difference of two quantum 
harmonic oscillators respectively in Z_ and Z + , where the role of h is taken by our 
constant A. The eigenstates of TC : 



Hil;{Z + ,Z-) = E^(Z + ,Z 
can be easily obtained. They are ip(Z + , ZJ) = ip° sc (Z + )ip. 



Z C {Z.) with 



MZ±) = (^ n n\o Q )~ l/2 H n [ — )exp 



zl 

2a 2 ) ' 



n 



0,+l,+2, 



(D.ll) 



(D.12) 



where H n are the Hermite polynomials and oq = J—- The eigenvalues are: 



Em 



A 



m H — ^ Aoj 
2J 



n 



(m — n)o; = iVa; 



(D.13) 



where can take every positive or negative integer value: N = 0, ±1, ±2, • • •. This 
confirms the discretization phenomenon we found before. Let us notice that the quan- 
tity A disappears from the spectrum, so it is just an artifact of the (q, X p ) representation 
and it is needed in ( p.9| ) only for dimensional reasons. There was no need of it in the 
first derivation, eqs. ( p.3|) -( |LT7D , of the discretization phenomenon. Second let us no- 
tice that, due to the difference of the two oscillators quantum numbers m and n above, 
the zero-point "energy"^] is zero differently than in the quantum case. Note also that 
there is an oo-order degeneracy in the sense that associated to the eigenstate E = Noj 
there is the set of eigenf unctions: ip = ■ip n (Z + )ip n+ N(Z^), where n can be any integer 
if N > while n > -N if iV < 0. 

This doubling of oscillators in the classical case is basically due to the fact that the 
classical KvN wave functions depend on a number of variables (q, X p ) which is double 
with respect to the quantum case. 

11 We have put quotation marks around the word "energy" because, as we explained in section 4, 
the E is not the energy but one of the eigenvalues of the evolution operator. 
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E Appendix 



In section 4 we have analyzed the properties of the Landau problem in the Schrodinger 
picture of the KvN formalism. In this appendix we want to study the same problem 
in the corresponding Heisenberg picture^. In particular we want to find out which 
are the constants of motion, i.e. the operators that commute with the generator of 
the time evolution 7i. These operators will give us some indications concerning the 
trajectory of the classical particle in a constant magnetic field. 

Let us remember the form of the Liouvillian in the Landau problem: 



H = —X x p x H ( X y + -B\ Px )[py- —x ) H \ z p. 



m 



m 



m 



(E.l) 



Defining v y 
relation: 



m 



p y Ay ) and noticing that [x,H] = ip x /m, we get the following 



[Vy,H] 



m 



[Bx, X x )p x = - 7 ^ r p x 



m 2 c 



m 2 c 



(E.2) 



eB 



If we introduce the Larmor frequency: to = — we can then easily prove that 

mc 

xq = x + Vy/u is a constant of motion. In fact, using (|E.2|) we get: 

ieB 



[x ,H] 



x + ^,n 

to 



i mc 

—Px + —5 

m eB 



m 2 c 



p x =0 



(E.3) 



In the same way the commutators of y and v x with the Liouvillian are: 



m 



eB_ 

Py X 



[v x ,H] 



ieB 
m 2 c 



eB 

Py X 



and so we obtain that also y = y — v x /u commutes with 7i: 



(E.4) 



mc, 



[y 07 H] = [y,H]-—[v x ,H] = 



(E.5) 



Now classically a particle in a constant magnetic field directed along z describes an 
helicoidal orbit whose projection on the x, y-plane is a circumference with a radius 
equal to the Larmor one QL ar : 



I)' 



(E.6) 



Using eqs. ( [E.2| )-( [ET4]) it is possible to prove that also the Larmor radius is a constant 
of the motion: 

\QlarM = - 2 [vlH} + - 2 [v 2 r H] = (E.7) 



12 All the objects appearing in this appendix are operators. Therefore we will not use explicitly the 
hat-symbol (" A ") to indicate them. 
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The Larmor radius can be written also in terms of the x, x$ and y, y$ operators in the 
following way: 

Qlar = + v 2 y ) = \{uj{y - y )] 2 + \{u{x Q - x)f = (x - x Q ) 2 + (y - y ) 2 (E.8) 

UJ UJ UJ 

Therefore (xo, yo) is the center of a circumference which is the projection of the orbit 
of the particle onto the plane (x,y) and g Lar is the corresponding radius. Note that 
in the KvN operatorial formalism the operators xq and yo are suitable combinations 
of the if operators and they commute among themselves. This implies that they can 
be determined with arbitrary precision. In quantum mechanics, instead, one can prove 
that the following relation holds: 

[x ,y ] = ^ (E.9) 
en 

and therefore, differently than in classical mechanics, there is an uncertainty relation 
involving the coordinates of the center of the circumference. 
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Figure Caption 

Figure 1: Aharonov-Bohm geometrical set up. 

Figure 2: Zeros of Bessel functions: m=l (continuous line), m=0.9 (dashed line). 



37 



Figure 1: Aharonov-Bohm geometrical set up 
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Figure 2: Zeros of Bessel functions: m=l (continuous line), m=0.9 (dashed line) 
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